Comment on "A generalized Langevin formalism of complete DNA melting transition' 
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We show that the calculated DNA denaturation curves for finite (Peyrard-Bishop-Dauxois (PBD) 
chains are intrinsically undefined. 



PACS numbers: 87.15.Aa,87.15.He,82.39.Pj 

In Ref. [H, the authors claim to solve a problem of 
equilibrium statistical physics by increasing damping in a 
Langevin simulation, i.e. by changing the dynamics. This 
is a fundamental flaw which though forms the basis of the 
work. Increasing P may change how fast the phase space 
is explored, but, if the simulations are performed prop- 
erly, it cannot alter the equilibrium distribution. The 
modification introduced by the authors apparently solves 
the problem of the divergence of (y) as the simulations 
start from a closed initial condition instead of using a 
set of initial conditions distributed over the whole phase 
space. The large increase of damping considerably slows 
down the exploration of phase space in the calculation, so 
much that the system stays in the vicinity of the closed 
initial configuration. 

Moreover, the dynamical argument used by the au- 
thors does not stand examination because it underesti- 
mates the actual time-scale which is provided by the ex- 
periments. Using the standard diffusion equation (with 
corrected prefactor compared to [H) 

P{y,t;yo) = {^^^^y^\M-mT{y~yof/ikBTt), 

(1) 

the authors state that "the probability of attaining in- 
finitely large values of the bead displacements is effectively 
zero". Although strictly speaking correct, it is easy to 
show that this does not prevent divergence on the exper- 
imental timescale, which is typically 2 min. per temper- 
ature step. Using the variables P = 0.5 ps~^, m = 350 u, 
and T = 300 K, the chance that the two beads are within 
reasonable distance, e.g < 100 A , is still large when the 
timescales of computer simulations are considered, i.e. 
« 94% for t = 1 ns. However, taking the experimentally 
relevant scale of t = 2 min. , this probability has dropped 
to less than 2 • 10"*^. 

A Monte Carlo study of the same systemj^] also suffers 
from this problem as we showed in a recent paper The 
large quantitative differences between [l| and should 
have alarmed the authors that at least one of these stud- 
ies reports unconverged results. Instead, the results re- 
ported in P] are found to be in remarkable agreement 
with the former study and experimental results [3]. 
In principle, both Langevin [H and MC [a] should un- 
surprisingly sample the same distribution. Yet, although 
the curves look qualitatively similar, for some tempera- 



tures they show values that differ by more than 40%. As 
a general trend, when compared to experiments Q the 
curves are shifted towards lower values along the tem- 
perature axis. A recent comment 

[1 on 0, showed that 
denaturation curves, calculated in this way, continue to 
shift to lower temperatures when the simulation time is 
extended. The apparent signature of premelting in [li@| 
that has been observed experimentally Q, is actually a 
numerical artifact in the theoretical studies. The present 
PBD model is unable to reproduce this experimental fea- 
ture as can be concluded from more accurate numerical 
studies d, [1]. 

The misconceptions appearing in this work and others 
on the same subject show that it is necessary to clar- 
ify some essential principles regarding the modeling of 
thermal DNA denaturation and on the computational 
methods used to address the problem. The Peyrard- 
Bishop-Dauxois (PBD)[^ model has some fundamental 
limitations which seem to be ignored It was not de- 
signed to sample the phase space corresponding to fully 
melted DNA molecules. Therefore, computing exper- 
imental quantities such as p (fraction of fully melted 
molecules) does not make sense for several reasons. 

First, as discussed above, the phase space that com- 
prises the fully denaturated molecules is infinitely larger 
than that of the (partially) closed molecules. As result, 
the DNA molecule will always visit and reside in the fully 
denaturated state at any temperature. The reason that 
this does not appear in the experiments is because there 
is an equilibrium between melting and self assembly in 
DNA solutions. Conversely, the PBD model describes a 
single DNA molecule in an infinitely diluted system. 

Secondly, each base is linked to its complementary base 
by a given on-site potential, determined by the initial 
state. Therefore assembly /disassembly reactions are un- 
physically favored in a particular pathway. Equilibrium 
intermediate states are not well sampled as mismatches 
are not considered, while these structural events play a 
very important role in the premelting phase, specially in 
bubble nucleation and dynamics. 

In spite of its limitations the model can give useful 
results in the temperature range where the molecule is 
not fully denaturated (dsDNA ensemble [31). Surpris- 
ingly, the denaturation of infinite DNA chains can also 
be properly described. The presence of this other infin- 



2 



ity, besides that of the unbounded phase space, ensures 
that the molecule remains in the double stranded state 
at all times This can be understood via the follow- 
ing reasoning: when moving in one direction along the 
chain, the sequence y^, yi+i, . . . basically represents an in- 
finite onc-dimcnsional random walk which must always 
visit the origin (assigning the whole molecule as being 
'closed') according to Polya's theorem. 

Moreover, recent technological progress ^ makes it 
possible to isolate the experimental signal of partially 
closed molecules from the fully denaturated popula- 
tion. For this reason we introduced the bias potential 
Vbias(ymin) which makcs it possible to study the PBD 
model at elevated temperatures. In particular, we fo- 
cused on the denaturation curves which describe the 
fraction of open base-pairs as function of temperature 
within the population of closed molecules. This property 
is well defined and gives a more direct comparison with 
experiments [4J than the phenomenological equation tak- 
ing into account the fully open molecules Q, which re- 
quires additional fitting parameters. They are hard to 
determine properly because they appear in exponential 
terms and the results are very sensitive to their values. 

The criticism on this bias potential Vbias(ymin) is un- 
founded. The article claims that there is no rational jus- 
tification of this particular form, but, in fact, Vbias can be 
any kind of steep increasing function. The exact shape 
has no influence on the results as it only affects the tail 
of the distribution. Our results were actually not car- 



ried out by MD or MC, but using a much more efficient 
transfer integral approach |8| for which Vbias is simply 
an infinite step-function. Also the second argument that 
it lowers the equilibrium values of / ( fraction of open 
base pairs considering all molecules), p and I is untrue. 
/ and p are ill-defined anyway and therefore we did not 
report on them. However, 1{T) is perfectly defined and 
remains exact when the bias potential is invoked. As 
Vbias (2/min) Only cxcludcs parts of phase space that does 
not contribute to Z, it has no effect on it. This can easily 
be verified mathematically but has also been confirmed 
by unbiased MC simulations 

To conclude, the simulation approach applied in this 
work [l[ and other related studies, e.g. is not well 
founded in terms of statistical physics and uses the PBD 
model in situations where it can not be expected to be 
valid. The apparent signatures of premelting in P, Q 
are nothing more than a memory effect of chosen initial 
conditions and do not hold when they are scrutinized by 
more accurate calculations [1, Q . 
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